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Using a recently introduced method to quantify the time varying lead-lag dependencies between 
' pairs of economic time series (the thermal optimal path method), we test two fundamental tenets 

of the theory of fixed income: (i) the stock market variations and the yield changes should be anti- 
, correlated; (ii) the change in central bank rates, as a proxy of the monetary policy of the central 

' bank, should be a predictor of the future stock market direction. Using both monthly and weekly 

I data, we found very similar lead-lag dependence between the S&P500 stock market index and the 

yields of bonds inside two groups: bond yields of short-term maturities (Federal funds rate (FFR), 
3M, 6M, lY, 2Y, and 3Y) and bond yields of long-term maturities (5Y, 7Y, lOY, and 20Y). In all 
cases, we observe the opposite of (i) and (ii). First, the stock market and yields move in the same 
direction. Second, the stock market leads the yields, including and especially the FFR. Moreover, 
■ we find that the short-term yields in the first group lead the long-term yields in the second group 

' before the financial crisis that started mid-2007 and the inverse relationship holds afterwards. These 

• , results suggest that the Federal Reserve is increasingly mindful of the stock market behavior, seen 

Ch . at key to the recovery and health of the economy. Long-term investors seem also to have been more 

C^_^ ' reactive and mindful of the signals provided by the financial stock markets than the Federal Reserve 

itself after the start of the financial crisis. The lead of the S&P500 stock market index over the 
bond yields of all maturities is confirmed by the traditional lagged cross-correlation analysis. 

PACS numbers: 89.65.Gh, 89.75.Fb, 89.75.Kd 
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^ : I. INTRODUCTION 

. Financial markets play a more and more important role in the economic system. Many financial variables have 
predictive power for output or inflation of the real economy. Financial markets are becoming increasingly important 
to the real economy due to their impact on output growth and inflation, among others [H-Q. As an important part of 
financial markets, stock markets can be considered as economy barometers 0, Q- As a consequence, monetary policy, 
which is usually based on inflation target and sometimes unemployment goals, is not independent of stock markets. 
There is a large financial economic literature concerned with the impact of and relationship between the monetary 
] policy of central banks and the performance of stock markets. The common wisdom asserts that (i) the stock market 
■ ■ variations and bond yield changes should be anti-correlated and (ii) the change in short-term interest rates, as a 
proxy of the monetary policy of the central bank, should be a predictor of the future stock market direction. The first 
assertion reflects the impact of capital cost on economic growth. The second statement is a corollary of the causal 
effect of the former one. 

Some of the most relevant results for our study that were obtained by previous scholars on these two statements 
include the following. Tobin's portfolio selection theory explained the stock price increases observed in times when 
the interest rate goes down as due to investors' preference for the higher yield of stock markets. Rigobon and Sack 
[lo| documented that an increase in short-term interest rate results in a decline in stock prices and in an upward shift 
and flatter yield curve. Bernanke and Kuttner [Tl1 | found that a hypothetical unanticipated 25-basis-point cut in the 
FFR target is associated with approximately a 1% increase in the broad stock indexes. Bj0rnland and Laitemo (T2j 
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found a significant relationship, which is however the inverse of (i) and (ii): a one percent increase of the stock market 
leads on average to a 4-basis-point increase of the interest rate. Two of us have also previously found that the stock 
market seems to influence the FFR, during the 2000-2003 US stock market antibubble [iMi. 

Here, using an extension of the so-called TOP technique [l3l - [l6| for the joint analysis of pairs of time series, we 
revisit the pertinence of these two assertions (i) and (ii) by estimating the lead-lag structure between the US stock 
market proxied by the S&P 500 index and a set of Treasury bond yields, including the Federal funds rate (FFR), which 
constitutes one of the tools implementing monetary policy in the US. Our analysis is applied to monthly and weekly 
data of Federal funds effective rate (FFR), and nine Treasury bond yields with different maturities: 3M (3 months), 
6M, lY (1 year), 2Y, 3Y, 5Y, 7Y, lOY, and 20Y. The period of analysis from August 2000 to February 2010 includes 
the bearish market up to mid-2003, the bullish bubble-like market regime up to October 2007 followed by the turbulent 
phases associated with the so-called great Recession [3 • Given the extraordinary developments associated with 
the financial crises followed by economic crises in different parts of the world, it is particularly interesting to investigate 
the lead-lag structure between the US stock market proxied by the S&P 500 index and a set of Treasury bond yields. 

This work is organized as follows. In Section HH we brief review the thermal optimal path method and propose to 
assess the statistical significance of the extracted lead-lag structure by using bootstrapping techniques. Section IIIII 
describes the data sets and presents the results of unit root tests for these time scries. Our main results arc presented 
in Section IIVI Comparisons with the traditional lagged cross-correlation analysis are given in Section |Vl Section IVII 
summarizes our findings. 



II. METHODOLOGY 



A. Description of the thermal optimal path (TOP) method 

The thermal optimal path (TOP) method has been proposed as a new method to identify and quantify the time- 
varying lead-lag structure between two time series. The TOP method was successfully applied to several economic 
cases [14| - [l6| . It works as follows. 

Consider two standardized time serie s \ X(t^) : ti — 1,...,A^} and {^^(^2) : ^2 = 1,...A^}. The matrix Ex,y of 
distances between A to F is defined as [ij, [l^| 

e{h,h) = [X{h)-Y{h)]^ . (1) 

The element e(<i,t2) of the matrix Ex,y thus compares the realization A(<i) of X at time ti with the realization 
Y{t2) of Y at time t2- The value \X{ti) — Y{t2)Y defines the distance between the realizations of the first time series 
at time ti and the second time series at time t2- The N x N matrix Ex,y thus embodies all possible point- wise 
pairwise comparisons between the two time series. Note that the distance matrix could be modified to deal with two 
non-monotonic time series, for which the TOP algorithm is essentially the same [l6| . 

Once the matrix Ex,y with elements given by Eq. ([1]) is obtained, an optimal path is determined that quantifies 
the lead-lag dependence between the two time series. Figure [T] gives a schematic representation of how lead-lag paths 
are defined [ij]. The first (resp. second) time series is indexed by the time ti (resp. 12). The nodes of the plane 
carry the values of the distance for each pair (ti,t2)- The path along the diagonal corresponds to taking ti = t2, i.e., 
compares the two time series at the same time. Paths above (resp. below) the diagonal correspond to the second 
time series lagging behind (resp. leading) the first time series. The figure shows three arrows which define the three 
causal steps (time flows from the past to the future both for ti and t2) allowed in our construction of the lead-lag 
paths. A given path selects a contiguous set of nodes from the lower left to the upper right. The relevance or quality 
of a given path with respect to the detection of the lead-lag relationship between the two time series is quantified by 
the sum of the distances along its length, called the "cost" of the path. The lead-lag structure is then obtained as the 
relationship t2{ti) as a function of ti, as described shortly. We stress that the two-layer scheme presented in Fig. [1] 
performs better than multi-layer schemes [l^. 

As shown in Fig. [1] it is convenient to use the rotated coordinate system [x, t) such that 

ti = l + {t-x)/2 
t2 = l + {t + x)/2 ' 

where t is in the main diagonal direction of the {ti, 12) system and x is perpendicular to t. The origin [x = Q,t — Q) 
corresponds to (ti = l,t2 = 1). Then, the standard reference path is the diagonal of equation x = 0, and paths which 
have x(t) 7^ define varying lead-lag patterns. Inverting we have 



X ^t2-ti. 



(3) 
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FIG. 1: (Color online) Representation of the two-layer approach in the lattice (ii,i2) and of the rotated frame {t,x) as defined 
in the text. The three arrows depict the three moves that are allowed to reach any node in one step. 



This means that a positive x corresponds to <2 > ii, which by definition of the optimal thermal path below means 
that the second time series Y{t2) lags behind the first time series X{ti), or equivalently X{ti) leads Y{t2)- 

The idea of the TOP method is to identify the lead-lag relationship between two time series as the best path in a 
certain sense. A natural idea is that the best path is the one which has the minimum sum of its distances along its 
length (paths are constructed with equal lengths so as to be comparable) . This path with minimum cost has thus the 
minimum average distance between the two time series, i.e., it is such that Y{t2) resembles the most X{ti) along this 
path ^2(^1)- The problem with this idea is that the noises decorating the two time series introduce spurious patterns 
which may control the determination the path which minimizes the sum of distances, leading to incorrect inferred 
lead-lag relationships. It has been shown that a robust lead-lag path is obtained by defining an average over many 
paths, each weighted according to a Boltzmann-Gibbs factor, hence the name "thermal" optimal path method |14l4l6l |. 
Intuitively, this corresponds to performing an averaging operation over neighboring paths of almost the same cost. 

Concretely, we first calculate the partition functions G{x,t), for all values of a; at a fixed t in the lattice shown in 
Fig. [U and their sum G{t) = J^x G{x,t) so that G{x,t)/G{t) can be interpreted as the probability for a path to be 
at distance x from the diagonal for a distance t along the diagonal. This probability G{x,t)/G{t) is determined as a 
compromise between minimizing the mismatch or cost as defined above (similar to an "energy" ) and maximizing the 
combinatorial weight of the number of paths with similar mismatches in a neighborhood (similar to an "entropy"). 
As illustrated in Fig. [T] in order to arrive at (ti -I- 1,^2 + 1), a path can come from {ti + 1,^2) vertically, (^1,^2 + 1) 
horizontally, or (^1,^2) diagonally. The recursive equation on G{x,t) is therefore 

G{x, t + l) = [G{x -l,t) + G{x + 1, t) + G{x, t - l)]e-'(^'*)/'^, (4) 

where e{x,t) is defined by Eq. ([T]). The parameter T plays the role of a "temperature" controlling the relative 
importance of cost versus combinatorial entropy. The larger T is, the larger the number of paths that contribute to 
the partition functions. In contrast, as T 0, only the path with minimum cost counts. The recursion relation 1^ 
is derived following the work of Wang et al. To get G(x, t) at the t-th layer, we need to know and bookkeep the 
previous two layers from G{-,t — 2) to G{-,t — 1). After G{-,t) is determined, these values are normalized by G{t) so 
that G{x, t) does not diverge at large t. The boundary condition of G{x, t) plays an crucial role. For t = Q and t = 1, 
G{x,t) = 1. For t > 1, the boundary condition is taken to be G{x = ±t,t) ^ 0, in order to prevent paths to remain 
on the boundaries. 

Once the partition functions G(x,i)'s have been calculated, we can obtain any statistical average related to the 
positions of the paths weighted by the set of G(a;,t)'s. For instance, the local time lag {x{t)) at time t is given by 

(x(t))==5]x-G(x-,t)/G(i) . (5) 

X 

Expression ([5]) defines {x{t)) as the thermal average of the local time lag at t over all possible lead-lag configurations 
suitably weighted according to the exponential of minus the measure e(a;, t) of the similarities of two time series. For 
a given xq and temperature T, we determine the thermal optimal path {x{t)). We can also define an "energy" or cost 
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eT(a^o) to this path, defined as the thermal average of the measure e{x,t) of the similarities of two time series: 

2Ar-l-|a:o 



B. Bootstrapping tests and statistical significance 

In order to test whether the extracted lead-lag structure is statistically significant, we introduce a bootstrap approach 
j20j that is specifically adapted to the present problem. This statistical test extends and make more robust the method 
and results, as compared with previous works [l3 - [T6| . Consider two time series X{ti) (for instance the logarithmic 
returns of S&P 500) and Y{t2) (for instance the time increments of bond yields). We perform the TOP analysis on a 
fixed time interval at some temperature T. Let us assume we obtain the lead-lag function x{t). Recall that t is the 
diagonal of the ti x t2 plane. We then shuffle X{t) and Y{t) and redo the TOP analysis at the same temperature 
T. We obtain a new lead- lag function xi{t). This process is repeated another n ~ I times, giving a total of n paths 
Xi{t) with i ~ 1,2, . . . ,n. A typical value of n used below is 1000. For each t, out of the n ~ 1000 reshuffled time 
series, we determine the 5% quantile 3:51% (i) and the 95% quantile Xg5%(t), denoted in the following as L(t) and U(t). 
If x{t) is smaller than L{t) or larger than U{t), we interpret that the lead-lag x{t) at time t is different from zero 
at the significance level of 95% or larger. Complementarily, given the obtained lead-lag x{t), out of the n = 1000 
reshuffled time series, we obtain the p- value as a function of t, which thus characterizes the time periods when there 
is a statistically significant lead-lag structure as those with small p- values. 

III. DESCRIPTION OF THE DATA SETS 
A. Data sets 

In the following, we apply the TOP method respectively to monthly and weekly data of the S&P 500 index, Federal 
funds effective rate (FFR), and nine Treasury bond yields with different maturities: 3M (3 months), 6M, lY (1 year), 
2Y, 3Y, 5Y, 7Y, lOY, and 20Y. Each time series spans from August 2000 to February 2010. The Treasury bond at 
30-year maturity is not considered because it was discontinued in January 2002 and then reintroduced in February 
2006. 

Figure [2] shows the weekly sampling of the FFR, the nine Treasury bond yields with different maturities, and the 
S&P 500 index. In the left panel of Fig. [5J very interesting patterns emerge in the term structure. In general, the 
yields of Treasury bonds with short maturities are more sensitive to the economic circumstance and change to a larger 
extent. In 2000, 2006 and 2007, the spread is very narrow and the FFR is even higher than the Treasury bond yields, 
corresponding to an anomalous inverted yield curve. These time periods correspond respectively to the early stages 
of the 2000 US stock market crash and to the current financial crisis. The spread reaches local maxima in 2004 and 
2010. In addition, the right panel of Fig. [5] suggests that the FFR and the S&P 500 index change roughly in the 
same direction. It is thus interesting to refine this visual impression and determine rigorously using the TOP method 
described above what is the lead-lag structure between the evolution of the US stock market and the FFR, which 
embodies an important part of the policy of the Federal Reserve. 

In this paper, we use as inputs the logarithmic returns of the S&P 500 index and the increments of the FFR and of 
all the yields, rather than the non-stationary original time series. We define the logarithmic returns of the S&P 500 
index as follows 

X{t) = \^\ShP(t)\ - \^\ShP(t - 1)], (7) 

and the logarithmic increments of yields curves as follows 

Y(t) = ln[yield(t)] - ln[yield(t - 1)] , (8) 

where the time unit for t is one week for weekly data and one month for monthly data. We then normalize the two 
time series X{£) and Y(t) so that their mean is zero and their standard deviation is equal to 1 [l3|- This ensures that 
they are comparable and can be used meaningfully in the TOP analysis to extract their lead-lag structure. 
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FIG. 2: (Left) Weekly sampling of the Federal effective funds rate (FFR) and nine Treasury bond yields. (Right) S&P 500 and 
FFR together with the 20Y for comparison. 



B. Unit root tests 



We perform unit root tests on the logarithms of the original time series and their first-order differences {X{t) 
and Y{t)) to check for their stationary. The augmented Dickey- Fuller (ADF) [2l[, Phillips-Perron (PP) [22|, and 
Kwiatkowski-Phillips-Schmidt-Shin (KPSS) [l^] tests are adopted. For the ADF and PP tests, the null hypothesis is 
that the time series has a unit root, which utilizes the t-statistic. In contrast, the null hypothesis of the KPSS method 
is that the time series is stationary and uses the LM-statistic. The results are presented in Table H 



TABLE I: Unit root tests of the logarithmic monthly and weekly data and their first-order differences. The augmented Dickey- 
Fuller (ADF), Phillips- Perron (PP) and Kwiatkowski-Phillips-Schmidt-Shin (KPSS) tests are adopted. is the critical value 
at the xYo significance level, V is the statistic, and p is the p-value. 





K% S&P 500 


FFR 


lY 


5Y 


20Y 


method 


Vio% V5% Vi% V p 


V 


P 


V 


P 


V 


P 


V 


P 


Logarithmic monthly data 


ADF 


-2.58 -2.89 -3.49 -2.14 0.23 


-0.56 


0.87 


-0.64 


0.86 


-2.27 


0.16 


-2.51 


0.12 


PP 


-2.58 -2.89 -3.49 -2.15 0.23 


-0.01 


0.95 


-0.23 


0.93 


-2.00 


0.27 


-2.51 


0.12 


KPSS 


0.35 0.46 0.74 0.16 >0.1 


0.31 


»0.1 


0.27 


>0.1 


0.43 


0.06 


1.02 


0.00 


Logarithmic weekly data 


ADF 


-2.57 -2.86 -3.44 -2.08 0.25 


-0.16 


0.94 


-0.29 


0.92 


-2.10 


0.21 


-2.95 


0.04 


PP 


-2.57 -2.86 -3.44 -2.08 0.25 


-0.00 


0.96 


-0.30 


0.92 


-2.04 


0.27 


-2.72 


0.07 


KPSS 


0.35 0.46 0.74 0.31 >0.1 


0.69 


0.01 


0.62 


0.02 


0.84 


0.00 


1.91 


0.00 


Difference of logarithmic monthly data 


ADF 


-2.58 -2.89 -3.49 -8.29 0.00 


-4.21 


0.00 


-6.68 


0.00 


-8.07 


0.00 


-9.35 


0.00 


PP 


-2.58 -2.89 -3.49 -8.35 0.00 


-5.36 


0.00 


-6.66 


0.00 


-8.07 


0.00 


-9.83 


0.00 


KPSS 


0.35 0.46 0.74 0.11 >0.1 


0.32 


»0.1 


0.32 


»0.1 


0.10 


>0.1 


0.05 


»0.1 


Difference of logarithmic weekly data 


ADF 


-2.57 -2.86 -3.44 -22.5 0.00 


-24.6 


0.00 


-16.6 


0.00 


-19.3 


0.00 


-17.4 


0.00 


PP 


-2.57 -2.86 -3.44 -22.5 0.00 


-24.6 


0.00 


-16.8 


0.00 


-19.4 


0.00 


-17.4 


0.00 


KPSS 


0.35 0.46 0.74 0.13 > 0.1 


0.35 


0.10 


0.31 


>0.1 


0.07 


>0.1 


0.04 


>0.1 



For the logarithmic monthly data and logarithmic weekly data, the ADF and PP tests show that these time series 
are not stationary and have a unit root since the p- values are greater than 10%, except for the 20Y yield. In contrast, 
the KPSS test suggests that four time series are stationary since the p-values are much greater than 10%. 

For the differences of the logarithmic monthly data and logarithmic weekly data, the ADF and PP tests show that 
all time series arc stationary at the 1% significance level, and the KPSS test also confirms that these time scries are 
stationary at the 10% level. These results justify our use of the logarithmic returns in the TOP analysis in order to 
avoid possible spurious signals in the estimated lead-lag structure that could result from large excursions exhibited 
by the non-stationary time series. 
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IV. THE S&P500 LEADS ALL YIELDS: EVIDENCE FROM THE TOP METHOD 

A. Empirical results 

Figure [3] shows the instantaneous evolution of the lead- lag x{t) between the returns of the S&P500 index taken as 
the first time series and the logarithmic variation of each of the yields for the monthly data at temperature T = 2. 
We have been careful to investigate the impact of the locations of the starting and ending extremities of the paths. 
There are indeed a total of 19 x 19 thermal optimal paths, because there are 19 starting points (ii, start, ^2, start) and 
19 ending points (ii.cnd, ^2,ond), denoted using the (^1,^2) system instead of the {t,x) system for simphcity. The 19 
starting points are {ti = 0, t2 = 0), (ti = 0, ^2 = and {ti = i, ^2 = 0) for i = 1, 2, . . . , 9. The 19 ending points are 
(<i = N,t2 = N), {ti = N,t2 ^ N - i) and {h ^ N - i,t2 = N) for i ^ 1,2, . . . ,9, where N is the length of the time 
series. The overall thermal optimal path x{t) is chosen as the one with minimal energy (or total cost) among the 
19 X 19 thermal paths. As for the choice of the temperature T, we investigated other values and found our results to 
be robust and qualitatively similar with respect to variations of T between 1 and 3. To present our results, we choose 
this value T = 2 as it seems to represent a reasonable optimal, confirmed by cross-correlation analyses performed on 
the steady periods found with fixed lag times for various T's. 




t t 

FIG. 3: (Color online) Dependence of the lead-lag x(t) between the returns of the S&P500 index taken as the first time series 
and the logarithmic variation of each of the yields for the monthly data: (a) FFR, 3M, 6M, lY, 2Y, and 3Y Treasury bond 
yields as the first group; (b) 5Y, 7Y, lOY, and 20Y bond yields as the second group. The unit of x{t) is one month. 

Figure [3] is organized in two panels, each panel plotting one group. The first group includes FFR, 3M, 6M, lY, 2Y, 
and 3Y Treasury bonds as shown in Fig. [3lja). The second group includes 5Y, 7Y, lOY, and 20Y Treasury bonds as 
shown in Fig. |3Jb). The evolution of x{t) in each group are quantitatively similar. 

Figure 0] is the same as Fig. [3] for weekly data. Apart from largest fiuctuations of the lead functions, the results 
are very similar and robust to this change of time scale from monthly to weekly. 

B. Statistical significance 

Before commenting and exploiting the information presented in Figs. [3] and 01 it is important to ascertain their 
statistical significance. For this, we use the bootstrap method described above in section III Bl Figure [5] illustrates 
the obtained results from the monthly data for two maturities, namely the shortest one (FFR) and the longest one 
(20-year Treasury bond yield). It shows that the two lead function x{t) are well above the 95% quantile curves, that 
is, x{t) > U{t). The conclusion is the same for other Treasury bond yields. We conclude that the obtained lead-lag 
structure for the monthly data cannot be produced by chance at the 95% significance level. 

Figure El illustrates the obtained results from the weekly data for two maturities, namely the shortest one (FFR) and 
the longest one (20- year Treasury bond yield). The conclusion is the same for other Treasury bond yields. Therefore, 
the x{i) functions for the weekly data are positive at the 95% significance level, which unveils the nontrivial intrinsic 
lead-lag structure of the S&P 500 index and the yield time series. 
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C. Two shocking stylized facts 

The first and most important observation extracted from Figs |3] and H] is that, for all yields and at all times, the 
S&P500 index leads the yield changes, since x{t) is always positive, which by definition ([3]), means that t2 > ti for 
the optimal thermal path. Since the index ti corresponds to the S&P500 index and the index t2 corresponds to one 
of the yields, this conclusion follows. This result confirms and extends considerably that reported previously by two 
of us [ij] using standard measures of correlations over a restricted period from 2001 to 2003, under the somewhat 
provocative title "Causal slaving of the U.S. Treasury Bond Yield ... by the Stock Market..." Indeed, as this title 
suggests, this result x(t) > is particularly striking and rich of implication. This result collides against the common 
wisdom that usually asserts the following two rules: 

(i) the stock market variations and the yield changes should be anti-correlated; 

(ii) the change in FFR, as a proxy of the monetary policy of the central bank, should be a predictor of the future 
stock market direction. 

Indeed, according to the standard story, a lower interest rate means lower costs of borrowing for the private sector, 
implying that the private sector is going to profit from this opportunity by increased investments in innovations and 
entrepreneurial opportunities, leading (with some lag) to an improved outlook for the future growth of the economy. 
Since stock market prices reflect the anticipation of investors, this better outlook for the future economy should be 
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t t 

FIG. 6: (Color online) Bootstrap test for the significance of the lead-lag structure for the FFR (a) and the 20Y Treasury bond 
yield (b). 

soon reflected in the appreeiation of the stock market. Reciprocally, an increase of the yields and in particular of the 
FFR should, according to the standard story, translate soon into a drag on the growth of stock markets. 

We observe the opposite of (i) and (ii). First, we find that the stock market and yields move in the same direction, 
as pointed out independently by R. Werner [2J|. Second, the stock market leads the yields, including and especially 
the FFR. The implication is clear: the central bank policy is (1) reacting to the stock market and (2) is following 
it. When the stock market exhibits a rally, the Fed tends to progressively increase its rates as an attempt to calm 
down the "overheating engine" , as occurred towards the end of the ICT bubble when the Fed rate was increased to 
6.5%. A similar increase of the Fed rate occurred from 2004 to 2007. When the stock market plunges, the Fed tends 
to decrease its rates, in the hope of putting a brake on the stock market losses that negatively feedback onto the real 
economy via the wealth effect. 

Both previous and present Fed chairmen Greenspan and Bernanke have increasingly made clear that the 
Federal Reserve does care more and more about the evolution of the stock markets. On Dec. 3rd, 2010, 
former Federal Reserve Chairman Alan Greenspan told CNBC that rising stock values have played a critical 
role in the economic recovery. The stock market got a boost from the Fed policy to boost liquidity, which 
drove interest rates down and pushed investors toward riskier investments like stocks. "I think we are un- 
derestimating and continuing to underestimate how important asset prices, very specifically equity prices, are 
not only to shareholders but the economy as a whole," he said. Equities have risen more than 80% from the 
lows set during the financial crisis, noted Greenspan, benefiting investors and helping fuel the recovery. [Source: 
http: //www.dailyf inance . com/ story/ investing/ greenspan-rising- stock-markets- are-key-to-recovery/ 19743325/?i 
On Nov. 3rd., 2010, Bernanke issued the following statement in an opinion article for the Washington Post released 
hours after the Fed announced the $600 billion of Treasury buying through June in a second round of unconventional 
monetary stimulus: "Resuming large-scale asset purchases should boost economic growth through lower borrowing 
costs and higher stock prices... Stock prices rose and long-term interest rates fell when investors began to anticipate 
this additional action... Easier financial conditions will promote economic growth." Being content to see the stock 
market growing, this suggests a hidden mandate of the Federal Reserve to steer the stock markets. 

It seems that the dynamics of the Fed policy, as translated in the Fed rates and the longer maturity yields (which 
of course are far from being controlled by the central bank), is much more straightforward than articulated in fancy 
models [2^ . The evidence presented here suggests that Fed policy appears to be as if a straightforward reaction to 
financial markets was the main factor. 



D. Comparison between different yields 

Comparing the lead functions x(t) for the various yields with different maturities, we find that the short-term yields 
in the first group (left panel of Figs. [3] and move approximately in synchrony with the long-term yields in the 
second group (right panel of Figs. [3]and|4]) until 2007. And this synchrony is almost perfect from the yields spanning 
FFR to 3Y in the first group until mid-2007. Thereafter, during the time period following the financial crisis that 
started in mid-2007, we can observe that the short-term yields clearly lead the long-term yields and we have the 
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sequence of inequalities 

2:FFR(i) > Xsuit) > Xeuit) > XiY{t) > X2Y{t) > Xsy^t) > . (9) 

This is seen from the fact that x(t) tends to be larger for the short-term yields, since they are all compared with 
the same S&P500 stock market index. It is also interesting to observe the increasing lag x{t) between the yield rates 
and the S&P500 index from around month in 2000 to about one year in 2007. This is followed by a plateau for all 
yields from FFR to 3Y, that last about 2 years and is then followed by a decay of the lag thereafter to about half its 
maximum, i.e., around 6 months. 

For the second group of yields with maturities from 5Y to 20Y whose x(i)'s are plotted in the right panel of Figs. 
|3]and|4j the picture is somewhat different. Before early 2003, the four curves are close to each other with no clear 
lead-lag structure between them. Then, from 2003 to mid-2007, a period corresponding to a very bullish upward 
trend of the stock market boosted by the favorable low rate of the Fed policy and a booming real-estate bubble, one 
can observe that the longer term yields lead clearly the shorter term yields: 

X5Y{t) < X7Y{t) < Xiorit) < X2GY{t). (10) 

Thereafter, in the reaction to the financial crisis, one observes as for the FFR-3Y yields that the shorter-term yields 
lead the long-term yields: 

X^Y{t) > X7Y{t) > Xiorit) > X20Y{t). (11) 

There is much less evidence for a plateau of the lead structure with respect to the S&P500. 

We would also like to mention that a reversal such as the one from ([TU)) to ([TT|) does not seem to have been 
documented before. 



V. THE S&P500 LEADS ALL YIELDS: EVIDENCE FROM CROSS-CORRELATION ANALYSIS 

By construction, the traditional cross-correlation analysis [1^ is not adapted to time varying lead-lag structures. 
It is however useful to investigate how it performs in the present context in which the TOP method has diagnosed a 
significant time varying structure. 

Two representative time series (FFR and 20Y) are presented for illustration. The significance levels of the cross- 
correlations are evaluated using bootstrapping tests through shuffling the return time series, similar to the analysis for 
the TOP method described in Section HlB. We use the monthly data in this analysis. For each pair of time series, we 
analyze the whole time series and two non-overlapping time periods. The results are shown in Fig. [71 It is obvious that 
the lagged cross-correlation analysis is not able to characterize the instantaneous evolution of the lead-lag structure 
evidenced in the previous TOP analysis. This is not a surprise. 



A. (S&P 500, FFR) pair 

For the (S&P 500, FFR) pair in the whole time period, the highest peak found in Fig. [7l[a), with a positive lag, 
shows that the FFR lags behind the S&P 500 index by about 3 months, with a cross-correlation coefficient c = 0.43, 
which is significantly positive at the confidence level of 95%. There are two other peaks that are also significant, one 
at a negative lag of i = — 2 month with c = 0.26 and another at the positive lag t = 8 month with c = 0.24. The 
largest peak with positive lag and highest cross-correlation coefficient c = 0.43 can be considered as confirming the 
main results of the previous section that the stock market changes precede the FFR variations. Due to the fixed 
lead-lag structure of the method, the cross-correlation provides only an average coarse representation of the real much 
richer and dynamical nature of the lead-lag structure. 

For the time period before April 2007, we see many peaks at positive and negative lags t that are significantly 
different from zero, as shown in Fig. El^b). It is hard to extract from this plot a clear picture about the lead- lag 
structure between the S&P 500 and FFR. In the presence of large variations of the lead-lag structure, it is not 
surprising that the cross-correlation analysis is not informative. 

For the time period after April 2007, we see in Fig. [Tj^c) a significant peak at the positive lag t = 8 month with 
c = 0.56. This lag is consistent in magnitude with the average value of the x{t) curve shown in Fig. [SJa). This clear 
signal in the cross-correlation analysis can be explained from the fact that the lead-lag has stabilized approximately 
above a value of 6 months, according to the analysis of the xit) function shown in Fig. [Hl^a) during the time period 
under investigation. 
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FIG. 7: Lagged cross-correlation between the logarithmic return of the S&P 500 index and the logarithmic difference of the 
FFR (a-c) and between the logarithmic return of the S&P 500 index and the logarithmic difference of the 20Y Treasury bond 
yield (d-f) during different time periods: (a,d) the whole time period from August 2000 to February 2010, (b,e) the time interval 
from August 2000 to April 2007, and (c,f) the time interval from May 2007 to February 2010. The ordinate axis shows the 
cross-correlation coefficients c{t). The unit of the lag time t along the abscissa is month. 



B. (S&P 500, 20Y) pair 

For the (S&P 500, 20Y) pair in the whole time period, there are two significant peaks around zero lag t = m 
Fig.EId). 

For the time period before April 2007, the signal is ambiguous although we can see several significant peaks in 
Fig. me). 

For the time period after April 2007, we see in Fig. ^i) only one significant peak at t = 3 month with c — 0.56. 
According to Fig. [5jb), the lead-lag x{t) decreases from about 6 to 1 month. Therefore, these two analyses give 
consistent results: on average, the S&P 500 index leads the 20Y Treasury bond yield by about 3 months. 

Comparing Fig. [7] for the cross-correlation analysis and Fig. [5] for the TOP analysis, we can conclude that the 
cross-correlation analysis can extract only part of the information and the TOP method is clearly superior. 

VI. CONCLUDING REMARKS 

In this work, we have adopted the thermal optimal path method to investigate the dynamics lead-lag structure 
between the S&P 500 index of the US stock market and Federal Funds rate, as well as several Treasury bond yields 
with different maturities. The time period that has been investigated runs from August 2000 to February 2010. Both 
monthly and weekly data have been used and we obtained consistent results. In all cases, the S&P 500 index is found 
to lead the FFR and the bond yields. This is quantified by the lead function x{t) found to be positive at a high 
statistical confidence level determined by bootstrapping tests. This finding is consistent with and extends significantly 
a previous work reporting that the US Federal Reserve was "slaved" to the stock market during the 2000-2003 US 
stock market antibubble |13| . 

According to the TOP analysis, we observed that the FFR and the Treasury bond yields can be divided into two 
groups. The first group contains FFR, 3M, 6M, lY, 2Y, and 3Y bond yields with short-term maturities and the 
second group contains 5Y, 7Y, lOY, and 20Y bond yields with long-term maturities. The lead functions x{t) between 
the S&P 500 index and the yields in each group have very similar quantitative shapes, while they are different at a 
quantitative level across the two groups. We found that the short-term yields in the first group lead the long-term 
yields in the second group before the current financial crisis around 2007 and the inverse relationship holds afterwards, 
namely the long-term yields lead the short-term yields after 2007. 
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For the first group, the lead function x{t) increases during tlie time period from 2000 to 2007, followed by a two- 
year-long plateau, and then plummets in late 2009. We also found that the yields (including FFR) with shorter 
maturity in the first group have a longer lag behind the S&P 500 index than for the longer maturities. In contrast, 
for the second group, the lead function x{t) increases till 2006 and then decreases. We observed a reversal of the 
order of the lead functions x{t) among the different maturities in 2007: a yield with shorter maturity has a shorter 
lag to the S&P 500 index before the reversal point and a yield with longer maturity has a shorter lag to the S&P 500 
index after the reversal point. Qualitatively, the reversal phenomenon is coincident with the outbreak of the current 
financial crisis. 

The lag of the FFR to the S&P 500 index can be interpreted in the light of comments of the previous and present 
Fed chairmen Greenspan and Bernanke that the growth of stock markets is "key" to the recovery and health of the 
economy. The evidence provided here suggests indeed that the FFR policy is in a significant part influenced by the 
recent past behavior of the stock markets (stock market — > Federal Funds rate). In plain words, the fact that the 
FFR follows the stock market direction can be interpreted as a direct attempt to limit its losses and revive it in times 
of bearish markets or to stabilize it in times of overly buoyant bubbling markets. 

As for the longer maturities, the lag structure with respect to the S&P 500 index reflects (i) a natural link in the 
term-structure that attach the longer maturities to the shortest maturity and (ii) the aggregate strategics of investors 
facing uncertainties over the long term behavior of the economy pTi. [28j. In the first sub-stage before 04/2007, we 
observe the causal relational flow from the stock market — > Federal Funds rate — >■ short-term yields — !■ long-term 
yields, and afterwards, we find the flow from the stock market — ?> long-term yields —> short-term yields — !> Federal 
funds rate. Thus, the lead-lag structures between the different yields changed after the financial crisis starting in 2007. 
This change can be rationalized by the strategies implemented by long-term investors in the face of growing global 
market uncertainties, such as central Banks of major Asian countries and pension funds which are heavily invested 
in the US long-term Treasury bonds [1^. The stern challenges faced by the US economy escalated the uncertainty 
which cascaded to exchange rate and inflation. Consequently, the long-term Treasury bonds became quite reactive 
to the behavior of stock markets, reflecting the actions of these long-term investors "flying to safety" : a plunge in 
the stock markets led to strong demand for the supposedly safe US Treasury bonds, pushing down mechanically the 
corresponding yields. This suggests that the long-term investors have been more reactive and mindful of the signals 
provided by the financial stock markets than the Federal Reserve itself after the start of the financial crisis. This may 
be due to the more complex agenda as well as the delicate role of the Federal Reserve, which has to take into account 
the impact of its interventions [25j . Caution and prudence on the part of the Fed in a time of high uncertainty may 
thus be the reason for this inversion of the lead-lag relationship between changes of yields of different maturities. 
However, the robust lead of the S&P 500 stock market index with respect to yields of all maturities remains the most 
important stylized fact unearthed by our study. 
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